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Abstract. The genus Trachylepis is widespread throughout most of continental Africa and its surrounding islands. How¬ 
ever, the majority of phylogenetic studies on this genus have focused on species occurring in eastern and southern Africa. 
Herein, we examine relationships among ten Trachylepis taxa that occur in Central and West Africa: T. affinis, T. albilabris, 
T. aureogularis, T. gonwouoi, T. maculilabris, T. mekuana, T. perrotetii, T. polytropis polytropis, T. polytropis paucisquamis, 
and T. quinquetaeniata. Five genes (two mitochondrial and three nuclear) were sequenced for 153 individuals, revealing 
much higher levels of diversity than previously realized, and justifying the need for future taxonomic investigations. Be¬ 
cause of high levels of morphological conservatism in Trachylepis, the taxonomy of each of these species is complex, and 
previously synonymized names may be available for several lineages. Molecular dating techniques suggest that while the 
two major clades of Trachylepis represented in this study diverged approximately 23 million years ago, the majority of di¬ 
versification took place in the last 17 million years. Further work is needed to fill in sampling gaps and increase genetic 
coverage for some clades before the full genetic diversity of this group can be realized. 

Key words. Squamata, Scincidae, phylogeography, taxonomy, species complex, Cameroon Volcanic Line, Albertine Rift, 
Kwahu Plateau, Pleistocene refugia, Miocene aridification. 


Introduction 

The scincid genus Trachylepis Fitzinger, 1843 is wide¬ 
spread throughout most of continental Africa and its sur¬ 
rounding islands and comprises over 80 species (Uetz et 
al. 2019). Trachylepis diversity in Central and West Afri¬ 
ca, with the exception of a few “hotspots,” including the 
Cameroon Volcanic Line in Central Africa (LeBreton 
1999, Mausfeld-Lafdhiya et al. 2004), has been consid¬ 
ered lower than elsewhere on the continent. However, the 
majority of phylogenetic research on these skinks has been 
conducted on East and South African or insular lineages 
(e.g., Lima et al. 2003, Gunther et al. 2005, Rocha et al. 
2010, PoRTiK et al. 2010, 2011, Sindaco et al. 2012, Vences 
et al. 2014), but limited studies on Trachylepis in Central 
Africa have hinted at higher levels of diversity than are cur¬ 
rently recognized (e.g. Mausfeld-Lafdhiya et al. 2004, 
Jesus et al. 2005, Chirio & LeBreton 2007, Ceriaco et al. 
2016, Marques et al. 2019, Weinell et al. 2019). 


Thirteen species of Trachylepis can be found across 
mainland Central and West Africa. We focus on ten taxa 
for which genetic sampling is available: T. affinis (Gray 
1838), T. albilabris (Hallowell, 1857), T. gonwouoi Al¬ 
len et al., 2017, T. maculilabris (Gray, 1845), T. mekuana 
(Ghirio & Ineich, 2000), T. perrotetii (Dumeril & Bib- 
RON, 1839), T. polytropis (Boulenger, 1903), T. quinque¬ 
taeniata (Lichtenstein, 1823), the subspecies T. polytropis 
paucisquamis (Hoogmoed, 1978), and the putative spe¬ 
cies T. aureogularis (Muller, 1885). Many of these species 
have widespread sub-Saharan ranges with the exception of 
T. mekuana, which is endemic to the Bamboutos Massif in 
the Cameroon Volcanic Line (Chirio & LeBreton 2007, 
Ceriaco et al. 2016). 

Most Trachylepis were referred to the genus Euprepes 
(or Euprepis) Wagler, 1830 during the 19* century, and 
to Mabuya (or Mabuia) Fitzinger, 1826 throughout the 
20* century. Mausfeld et al. (2002), in the context of the 
dismantling of circumtropical Mabuya, resurrected Eupre- 
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pis as the name applicable to the primarily African line¬ 
age. However, Bauer (2003) demonstrated that Trachylepis 
was the oldest available name for this clade. Because of 
high levels of morphological conservatism in this genus, 
many of these species have complex taxonomic histories. 
Trachylepis affinis was described by Gray in 1838 from a 
single museum specimen of unknown origin. The nomi¬ 
nal taxa Euprepes blandingii Hallowell, 1844, Euprepis 
raddoni Gray, 1845, Euprepis frenatus Hallowell, 1857, 
Euprepes aeneofuscus Peters, 1864, Euprepes gracilis Bo- 
cage, 1872, and Euprepes pantaenii Fischer, 1885 were all 
synonymized with T. affinis in the comprehensive taxo¬ 
nomic revision of Hoogmoed (1974), which simultane¬ 
ously raised Euprepis albilabris to species status. Trachy¬ 
lepis maculilabris was described from “W. Africa,” and has 
at various times been considered to include the subspecies 
Mabuia m. albotaeniata Boettger, 1913, Mabuya m. co- 
morensis Peters, 1854, Mabuya m. boulengeri Sternfeld, 
1911, Mabuya m. casuarinae Broadley, 1974, and Mabuya 
m. infralineata Boettger, 1913. Today all of these taxa are 
considered valid species (Uetz et al. 2019), whereas the one 
subspecies and six varieties of T. maculilabris from eastern 
Africa described by Sternfeld (1912) are now considered 
to be synonyms of the nominotypical form (Geri'aco et 
al. 2016). Trachylepis polytropis comprises two subspecies, 
the nominotypical form and T. p. paucisquamis, the latter 
taxon is sometimes considered as specifically distinct, e.g. 
by Trape et al. 2012, who also considered T. aureogularis as 
a distinct species, although many authors consider it as a 
synonym of T. albilabris (e.g., Uetz et al. 2019). 

Previous research has suggested that the widespread 
species Trachylepis maculilabris, T. affinis and T. perrotetii 
represent species complexes (Mausfeld-Lafdhiya et al. 
2004, Ghirio & LeBreton 2007, Geri'aco et al. 2016), but 
genetic analyses of these taxa are lacking, and T. polytropis, 
T. p. paucisquamis, T. aureogularis and T. mekuana have 
never been assessed in an explicit phylogenetic context. 
This study aims to conduct the most comprehensive phylo¬ 
genetic analysis of Central and West African Trachylepis to 
date, and to provide a comprehensive time calibrated as¬ 
sessment of historical biogeography for this group. 

Materials and methods 

Taxon sampling 

A total of 153 samples from the species Trachylepis affinis, 
T. albilabris, T. aureogularis, T. gonwouoi, T. maculilabris, 
T. mekuana, T. perrotetii, T. polytropis, T. p. paucisquamis, 
and T. quinquetaeniata were examined from through¬ 
out Central and West Africa. Two species from the genus 
Eutropis, the sister group to the rest of the circumtropi- 
cal Mabuyinae radiation (Karin et al. 2016), were used 
as outgroups. Specimens from Cameroon, Nigeria, Niger, 
the Gulf of Guinea islands, Ghana, Gabon, Guinea, Mali, 
Uganda, the Central African Republic, Democratic Re¬ 
public of the Congo, Burundi, Tanzania, and Mozambique 
were collected by the authors or borrowed from the fol¬ 


lowing institutions: American Museum of Natural History 
(AMNH), California Academy of Sciences (CAS), Monte 
L. Bean Life Science Museum (BYU), Museum of Com¬ 
parative Zoology, Harvard University (MCZ), Museum of 
Natural History of Geneva (MNHG), Museum of Verte¬ 
brate Zoology, Berkeley (MVZ), North Carolina Museum 
of Natural Sciences (NCSM), United States National Mu¬ 
seum of Natural History (USNM), University of Texas at El 
Paso Biodiversity Collections (UTEP), and the University 
of Washington Burke Museum (UWBM). A list of all spec¬ 
imens used in this study and associated GenBank numbers 
can be found in Supplementary Table Si. 

Laboratory protocols 

Five genes, two mitochondrial: 16S rRNA (16S) and NADH 
Dehydrogenase Subunit 2 (ND2); and three nuclear: Re¬ 
combination Activating Gene 1 (RAGi), Exophilin 5 
(EXPH5) and Kinesin Family Member 24 (KIF24), were 
sequenced for each specimen. Genomic DNA was isolated 
from liver tissue stored in 95% ethanol with salt extractions 
(Aljanabi & Martinez 1997). The target genes were am¬ 
plified using double-stranded polymerase chain reaction 
(PGR). A list of gene specific primers can be found in Ta¬ 
ble 1. All PGR reactions used an Eppendorf Mastercycler 
gradient thermocycler for 34 cycles with an annealing tem¬ 
perature of 50°C, followed by visualization of PGR prod¬ 
uct on a 1.5% agarose gel stained with ethidium bromide. 
PCRs were cleaned using a homemade magnetic bead so¬ 
lution (Rohland & Reich 2012). Gycle sequencing was 
performed using Big Dye V3.1, followed by an additional 
cleaning with the magnetic bead solution. Samples were 
sequenced using an ABI 3730 automated sequencer at Vil- 
lanova University. 

The 5’ and 3’ strands were sequenced separately, then 
combined in Geneious™ version 5.6 (Kearse et al. 2012), 
with ambiguous base pairs edited by eye. Sequences were 
aligned using MAFFT v.7 (Katoh et al. 2002) with de¬ 
fault parameters (gap opening penalty = 1.53, offset value = 
0.123). Uncorrected genetic distance matrices were gener¬ 
ated in MEGA 7 (Kumar et al. 2016) 


Gene tree analyses 

We concatenated mitochondrial and nuclear data and im¬ 
plemented PartitionFinder (Lanfear et al. 2014) to de¬ 
termine the appropriate substitution model for each gene 
using the greedy algorithm and the Bayesian Information 
Griterion (BIG) for model selection. Phylogenetic analy¬ 
ses were performed in a Maximum Likelihood framework 
using RAxML v.8.1.1 (Stamatakis 2014) with the GTR- 
GAMMA model. Support was assessed using bootstrap 
analysis (Felsenstein 1985) with 1000 pseudoreplications. 
We constructed a Bayesian phylogeny of our partitioned 
dataset with Mr. Bayes v. 3.2 (Ronquist et al. 2012). The 
Markov chain was run for 20 million generations for two 
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Table 1. Primers used to amplify the genes in this study. 


Gene 

Primer 

Sequence 

Citation 

ND2 

METFl 

5’-AAGCTTTCGGGCCCATACC -3’ 

Macey et al. 1997 

ND2 

COIRI 

5’-AGRGTGCCAATGTCTTTGTGRTT-3’ 

Arevalo et al. 1994 

16S 

16sA-L 

5’ CGCCTGTTTATCAAAAACAT-3’ 

Simon et al. 1994 

i6S 

16sB-H 

5’- CCGGTCTGAACTCAGATCACGT-3’ 

Simon et al. 1994 

RAGI 

RAGlSFf 

5’-TTCAAAGTGAGATCGCTTGAAA-3’ 

PoRTiK et al. 2010 

RAGI 

RAGiSR1200 

5’- CCCTTCTTCTTTCTCAGCAAAA-3’ 

PORTIK et al. 2010 

EXPH5 

EXPH5F1 

5’-AATAAACTKGCAGCTATGTACAAAACAAGTC-3’ 

PoRTiK et al. 2010 

EXPH5 

B1811R 

5’-CRCACGTCTAGAACCAAAGGTCC-3’ 

Karin et al. 2016 

KfF24 

KIF24F1 

5’-SAAACGTRTCTCCMAAACGCATCC-3’ 

PoRTiK et al. 2010 

KIF24 

KlF24Ri 

5’-WGGCTGCTGRAAYTGCTGGTG-3’ 

PoRTiK et al. 2010 


separate runs. Convergence of the posterior likelihood val¬ 
ues was assessed using Tracer v 1.6 (Rambaut et al. 2014). 
The first 25% of trees were discarded as burn in, and the 
phylogeny was visualized using FigTree 1.4.2 (Rambaut 
& Drummond 2012). Maps were created in QGIS (QGIS 
Development Team 2014) with Natural Earth: vector and 
raster map data at naturalearthdata.com. 

Species tree analysis 

We performed a species tree analysis and estimated diver¬ 
gence times between Trachylepis species using the program 
■•^BEAST implemented through BEAST v. 1.8.2 (Drum¬ 
mond & Rambaut 2007). Several taxa, including T. me- 
kuana, T. polytropis and T. p. paucisquamis, were exclud¬ 
ed from the species tree because of missing data for some 
genes (Supplementary Table Si). Heterozygous sites were 
phased using the program seqPHASE (Plot 2010) and 
PHASE 2.1.1 (Stephens et al. 2001, Stephens & Don¬ 
nelly 2003). Each gene was partitioned separately for a 
total of five partitions. The following substitution mod¬ 
els were applied to each partition based on BIG support: 
16S (GTR -iG), ND2 (TN93-1G), RAGi, EXPH5 and KIF24 
(HKY-i-G). For the species tree prior we specified a Yule 
model and for the population size we used a piecewise lin¬ 
ear model with a constant root prior. We used uncorrelat¬ 
ed lognormal relaxed clocks on each of the gene partitions 
and, as no fossils are available for Trachylepis or even the 
wider Mabuyinae group, divergence times were estimat¬ 
ed using both 16S and ND2 substitution rates. Variation 
in rates was specified by a normally distributed prior on 
the lognormal clock mean. 16S was given a mean rate of 
0.0080 with a standard deviation of 0.0020, and ND2 was 
given a mean rate of 0.00895 with a standard deviation of 
0.0025. These rates are based on Barley et al. (2015) and 
represent the span of estimated mitochondrial clock rates 
across many different lizard species. All other genes were 
assigned clock mean priors from o to 0.2 and standard de¬ 
viations with exponential distributions with a mean of 0.05 
(Barley et al. 2015). Substitution parameters were assigned 
a uniform prior from o to 100 (initial = 1) and the gamma 


shape parameters were assigned a uniform prior from o to 
10 (initial = 0.5). The Markov chain was run for 300 million 
generations with sampling every 30,000 generations. Con¬ 
vergence of the posterior likelihood values was assessed 
using Tracer v 1.6 (Drummond et al. 2012), with the first 
30 million trees discarded as burn in. The species tree was 
visualized in FigTree 1.4.2, and variation in the species tree 
was visualized using a cloudogram created in DensiTree 
V2.2.2 (BOUCKAERT & HeLED 2OI4). 

Results 

Phylogenetic analyses 

Our concatenated dataset consisted of 4376 base pairs: 16S 
(540 bp), ND2 and tRNA flanking region (1376 bp), RAGi 
(1095 bp), EXPH5 (786 bp), and KIF24 (579 bp). Our da¬ 
taset contained 26.9% missing data. PartitionFinder iden¬ 
tified the following models of substitution for each gene 
for the MrBayes analysis: 16S (SYM-i-I-i-G), ND2 cod¬ 
ing region (HKY-iI-iG) ND2 tRNAs (SYM-hI-hG), RAGi 
(HKY- iG), EXPH5 (HKY-iG), and KIF24 (K8o-tG). The 
model K80-1-G was not available in MrBayes and HKY-i-G 
was used instead. Topologies between the maximum likeli¬ 
hood and Bayesian phylogenies did not differ at any highly 
supported (>95% pp and 70 bs) node. Uncorrected genetic 
pairwise distances were 4.3-13.7 for 16S, 12.8-24.5 for ND2, 
1.3-7.3 for EXPH5, 0.4-5.5 for KIF24, and o.5-2.3 for RAGi 
(Supplementary Tables S2 and S3. These values are similar 
to those found for mitochondrial and nuclear loci in other 
skink genera (e.g. Poulakakis et al. 2005, Barley et al. 
2013, Rabosky et al. 2009). 

Trachylepis p. paucisquamis, which had previously been 
thought to be closely related to, or even a subspecies of, T. 
polytropis, was found to be part of a species complex that 
includes T. aureogularis and T. albilabris, whereas T. poly¬ 
tropis was found to be more closely related to T. maculi- 
labris (Fig. 1). Trachylepis affinis was recovered as a mono- 
phyletic group containing two divergent lineages, one from 
Central-West Africa and one south of the Kwahu Plateau 
in Ghana (Fig. 2). Trachylepis aureogularis was found to 
be distinct from, but sister to, T. albilabris. Trachylepis 


233 





Kaitlin E. Allen et al. 


r tb2.3 

—4 

I—fc 

I M 


-AMNH 109798 Trachytepis Quinquetaeniata 
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I MVZ 249759 Trachylepis quinquetaeniata 
I MVZ 249765 Trachyiepis quinquetaeniata 
p UWBM 6050 Trachylepis quinquetaeniata 
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Figure 1. Bayesian phylogeny of all Trachylepis species sampled dur¬ 
ing this study; well-supported nodes (PP > 0.95) are denoted by a 
black circle. 
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perrotetii was found to comprise three lineages: one from 
Guinea and northwestern Ghana, one widespread across 
Ghana, Niger and northern Cameroon, and one from the 
Cameroon Volcanic Line in northwest Cameroon and 
western Nigeria (Fig. 2). 

Trachylepis maculilabris was recovered as a species com¬ 
plex, more diverse than previously recognized (see Maus- 
feld-Lafdhiya et al. 2004), and comprising four previ¬ 
ously unidentified and divergent lineages. Of these lineag¬ 
es, one clade is distributed across Cameroon and Ghana, 
two are widespread throughout the Congo Basin, and one 
is isolated in the Albertine Rift in eastern Democratic Re¬ 
public of Congo (Fig. 2). Trachylepis mekuana, which has 
never before been placed in a phylogenetic context, was 
recovered as sister to the clade comprising T. perrotetii, 
T. gonwouoi, T. p. paucisquamis, T. albilabris, and T. affinis. 

Species tree 

The ’'^BEAST species tree analysis recovered similar rela¬ 
tionships to the maximum likelihood and Bayesian analy¬ 
ses, with the exception that it recovered Trachylepis gon¬ 
wouoi as sister to T. affinis instead of the clade containing 
T. albilabris, T. aureogularis, T. p. paucisquamis and T. affi¬ 
nis (Fig. 3). Uncertainty in the ’'^BEAST tree can be visual¬ 
ized with a DensiTree plot (Eig. 4). This plot shows high 
levels of uncertainty in the placement of T. gonwouoi and 


T. quinquetaeniata, and in the relationships between the 
relatively young T. maculilabris lineages. 

Dating analyses recovered the two major clades of 
Trachylepis in Central-West Africa diverging from one an¬ 
other during the late Oligocene, around 23.8 million years 
ago. Subsequent diversification in the T. gonwouoi, T. albi¬ 
labris, T. aureogularis, T. p. paucisquamis and T. affinis clade 
occurred during the mid- to late Miocene, and in the T. ma¬ 
culilabris clade during the Pliocene and Pleistocene (Fig. 3). 

Discussion 

Genetic diversity 

Several taxa of Trachylepis in Central and West Africa were 
included in this study that had never before been placed 
in a phylogenetic context. These species include T. aure¬ 
ogularis, T. mekuana, T. polytropis, and its subspecies T. p. 
paucisquamis. Trachylepis mekuana is a Cameroon en¬ 
demic found only on the Bamboutos Massif This species 
was found to be highly divergent and sister to the clade 
containing T. perrotetii, T. gonwouoi, T. p. paucisquamis, 
T. albilabris, and T. affinis. 

Based on morphology, Trape et al. (aora) suggested 
that T. aureogularis should be raised from synonymy with 
T. albilabris. Our results suggest that the two species are 
closely related but genetically distinct, providing further 
evidence in support of this interpretation. Likewise, Trape 



Figure 2. Sampling localities and major clades in A) Trachylepis perrotetii, B) the T. albilabris complex: T. albilabris (triangles), T. au¬ 
reogularis (circles) and T. p. paucisquamis (squares) C) T. maculilabris and D) T. affinis. The distributional ranges, shown in dark grey, 
are approximate and courtesy of the Global Assessment of Reptile Distribution Group. 
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et al. (2012) suggested that T. p. paucisquamis is morpho¬ 
logically distinct from T. polytropis and should be raised 
from its subspecies status. Our analyses support this assess¬ 
ment and found T. p. paucisquamis to be the sister species 
to a clade including T. albilabris and T. aureogularis, and 
T. polytropis is sister to T. maculilabris. While describing 
the subspecies T. p. paucisquamis, Hoogmoed (r974, r978) 
noted that it differed from T. polytropis in having a solid, 
versus broken, lateral band, supranasals in contact, and in 
the mid-body and sub-lamellar scale counts. Geographi¬ 
cally, these species are found in two distinct ranges: T. p. 
paucisquamis and T. aureogularis in West Africa, through¬ 
out the Ivory Coast, Liberia and Ghana; and T. polytropis 
and T. albilabris in Central Africa, in Cameroon, Bioko Is¬ 
land, Democratic Republic of Congo, Gabon, and the Cen¬ 
tral African Republic (Trape et al. 2or2, Uetz et al. 2or9; 
Fig. 2). 

Many of the species evaluated in this study contained 
several highly divergent lineages or comprised species 


complexes. An integrated taxonomic approach combining 
a more thorough sampling of genetic data, as well as mor¬ 
phology, is recommended before species diversity within 
these complexes can be fully assessed (Radial et al. 2010). 
However, high levels of unrecognized diversity would not 
be surprising given the morphologically conserved nature 
of Trachylepis and the paucity of research focusing on the 
Central African representatives of this group in the past 
(Mausfeld-Lafdhiya et al. 2004). 

Names may be available for several of these divergent 
lineages. For instance, Sternfeld (1912) described one 
subspecies and six varieties of Trachylepis maculilabris, 
most from the Albertine Rift of eastern Democratic Re¬ 
public of the Congo and Uganda. The varietal forms ber- 
geri, graueri, kwidjwiensis, rohrbecki, schubotzi, and wau- 
ensis, being infrasubspecific, are unavailable according to 
the International Code of Zoological Nomenclature (Bau¬ 
er et al. 2003, Ceri'aco et al. 2016), but the subspecies, 
T. m. major, considered a synonym of T. maculilabris, may 
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Figure 3. Dated *BEAST species tree including all major lineages of Trachylepis in this study. Divergence date ranges are shown in 
blue, and median divergence dates are listed by each node. Well-supported nodes (PP >0.95) are denoted by a black circle. 
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require reinvestigation. A divergent lineage of Trachylepis 
affinis south of the Kwahu Plateau in Ghana may corre¬ 
spond to one of several synonymized names from western 
Africa including: Euprepis raddoni, Euprepis aeneofuscus, 
Euprepis gracilis, Euprepis stangeri Muller, 1882, and Eu¬ 
prepis pantaenii. Additional divergent lineages that require 
further sampling as well as taxonomic assessment include 
the Trachylepis polytropis individual (UTEP 21830) from 
the Salonga River in Democratic Republic of the Congo, 
the T. gonwouoi individual from Ghana (MVZ 25261) and 
the T. perrotetii clade from Cameroon (Figs 1 and 2). 

Biogeography 

Although sampling areas were strongly biased (i.e., Gha¬ 
na often being the only sampled country in West Africa, 
no sampling in Nigeria and sampling only for Trachylepis 
maculilabris and T. polytropis in Democratic Republic of 
the Congo, Fig. 2), several common biogeographic patterns 
can be seen across the Trachylepis species in this study. The 


most common of which was in Ghana, where high levels 
of genetic divergence existed either across the Volta River 
(now Lake Volta), or south of the Kwahu Plateau in south¬ 
western Ghana (Fig. 2). This area was likely a rainforest 
refugium throughout the glacial cycles of the Pleistocene 
(Maley 1996), but the diversification seen between T. albi- 
labris and T. aureogularis is older, dating to the late Mi¬ 
ocene (Fig. 3). It is possible that the same processes that 
created rainforest refugia during the colder, drier glacial 
periods of the Pleistocene were also in effect during the 
cooling and aridification of sub-Saharan Africa during the 
Miocene (Maley 1996, Haffer 1997). 

High levels of endemism can be found within the Gam- 
eroon Volcanic Line. Two montane endemics, Trachylepis 
nganghae and T. mekuana were described from this area 
(Ghirio & LeBreton 2007). This study identifies another 
divergent lineage in the T. perrotetii complex occurring in 
an isolated region of the Gameroon Volcanic Line (Fig. 2). 
The T. perrotetii clade originated in the mid-Miocene, a 
period of dramatic cooling and major vegetation shifts 
throughout sub-Saharan Africa (Jacobs 2004, Kissling et 
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Figure 4. Cloudogram of uncertainty in the Trachylepis species tree. Individual trees (n = 5770) are shown in green and the consensus 
tree is shown in black. 
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al. 2012, Menegon et al. 2014), and may have diverged as 
a result of the uplift of the Cameroon Volcanic Line (Co- 
NARCCHiA & Dars 1983, Ubangoh et al. 1998, Bate Ti- 
BANG et al. 2017). 

Trachylepis maculilabris was the only species for which 
extensive sampling was available across the Democratic 
Republic of the Congo. Three clades were found in this 
country, two widespread throughout lowland areas and 
one in the Albertine Rift in the east. All of these diver¬ 
gences are relatively young and date back to the Pliocene 
and Pleistocene. It is possible that separate historic rain¬ 
forest refugia in Cameroon and the Congo played a role 
in the diversiftcation of these lineages (Maley 1996). Fur¬ 
ther analysis is necessary to determine if these divergences 
represent distinct species or genetically structured popu¬ 
lations. The Albertine Rift highlands house an extremely 
high number of endemic species (Greenbaum 2017), espe¬ 
cially birds (Stuart et al. 1990), amphibians (Greenbaum 
& Kusamba 2012) and reptiles (Hughes et al. 2018), and it 
is likely that some of the genetic diversity recovered here is 
representative of unrecognized taxonomic diversity. 

High levels of uncertainty were seen in the species tree 
for the placement of Trachylepis gonwouoi relative to T. af- 
finis, T. albilabris and T. aureogularis, and in the relation¬ 
ships between the T. maculilabris populations from Dem¬ 
ocratic Republic of Gongo (Fig. 4). In both cases, clades 
are either largely overlapping in geographic distribution or 
are parapatric (i.e., the Albertine Rift T. maculilabris). It 
is possible that historic introgression is responsible for the 
uncertainty in the relatively young clades (T. affinis, T. gon¬ 
wouoi, T. albilabris clade: 12 MYA; T maculilabris clade: 1.8 
MYA). It is also possible, particularly with T. gonwouoi and 
T. quinquetaeniata, that there was simply not enough in¬ 
formation in our matrix to accurately place these taxa, and 
further geographic and taxonomic sampling is necessary. 

Divergence dates 

The two major clades of Trachylepis in Gentral Africa di¬ 
verged from each other approximately 23.8 million years 
ago, around the end of the Oligocene. However, most of 
the diversiftcation within these clades occurred during the 
last 17 million years, from the mid-Miocene to the present. 
The period from the Miocene through the Pleistocene was 
characterized by dramatic climatic changes, including cy¬ 
cles of aridiftcation and humidiftcation across Africa (Ma¬ 
tey 1996, Matey & Brenac 1998) that correspond to di¬ 
versiftcation patterns in many other organisms (e.g., trees: 
Hardy et al. 2013; rodents: Nicotas et al. 2011, 2012; pri¬ 
mates: Gtifford et al. 2004, Tosi 2008, Haus et al. 2013). 
Several of the divergences between clades are distributed in 
areas that have been previously hypothesized as Pleistocene 
rainforest refugia (Matey 1996). These lineages include the 
three clades of T. maculilabris in Democratic Republic of 
the Gongo, and T. affinis and T. gonwouoi south of the Kwa- 
hu Plateau, suggesting that these areas of climatic stabil¬ 
ity may have played a role in rainforest species dynamics 


over a relatively long timeframe in Gentral Africa (Haffer 
1997). These results are also congruent with recent herpe- 
tological studies that suggest allopatric speciation associat¬ 
ed with rainforest refugia in the Upper and Lower Guinean 
forest blocks (Portik et al. 2017, Kpan et al. 2018), across 
Central Africa (Jongsma et al. 2018) and even in the north¬ 
west African Sahara-Sahel region (Gon^atves et al. 2018). 

This study highlights the high degree of unrecognized 
genetic diversity in Central African Trachylepis, but more 
work is necessary to fully understand the origins and ex¬ 
tent of this diversity and additional taxonomic implica¬ 
tions. Our data support Trape et al. (2012) in elevating 
T. aureogularis from synonymy and T. paucisquamis from 
subspecies status, and we recognize both taxa as valid spe¬ 
cies. Further taxonomic work is needed to determine if the 
additional genetically divergent lineages recovered here 
are morphologically diagnosable, and if any warrant the 
assignment of names currently residing in synonymy. As 
several of these lineages represent range-restricted endem¬ 
ics, conservation measures may need to be taken to assess 
the viability of populations and to protect them from the 
widespread deforestation in Central African rainforests. 
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